Masses and couplings of vector mesons from the pion electromagnetic, 
weak, and transition form factors 
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We analyse the pion electromagnetic, charged-current, and 717 transition form factors at timelike 
momentum transfers q, ~ s < 1.4 GeV^, using a dispersion approach. We discuss in detail the 
propagator matrix of the photon-vector-meson system and define certain reduced amplitudes, or 
vertex functions, describing the coupling of this system to final states. We then apply the derived 
analytic expressions to the analysis of the recent e^e~ — > tt^tt", t~ -K~-K^Vr, and e^e~ — > ■k'^^ 
data. We find the reduced amplitudes for the coupling of the photon and vector mesons to two 
pseudoscalars to be constant, independent of s, in the range considered, indicating a "freezing" 
of the amplitudes for s < 1 GeV. The fit to the form factor data leads to the following values 



of the Breit-Wigner resonance masses m 



775.3 ± 0.8 MeV, 



773.7 ± 0.6 MeV and 



782.43 ± 0.05 MeV, where the errors are only statistical. 



PACS numbers: 12.40.Vv, 13.40.Gp,13.65.-(-i 



I. INTRODUCTION 

The pion elastic and transition form factors at timelike momentum transfers provide an important source of in- 
formation about the masses and coupling constants of vector meson resonances p, u>, 0, etc. However a reliable 
extraction of the resonance parameters from the experimental data is a comphcated problem. The reason is that the 
direct QCD-based calculation of the form factors in terms of the resonance parameters is not possible yet, and hence 
one has to use approximate approaches. A typical procedure of extracting the vector meson parameters is as follows: 
One relies on some theoretical formula for the form factor in terms of the vector-meson masses and couplings, and 
tries to extract their numerical values by fitting the experimental data. Most of the approaches providing theoretical 
inputs for the form factors in the region of vector meson resonances 0.7 to 1.2 GeV fall into the two big classes: 
approaches based on vector meson dominance picture Q Q, 01 Q IH HI approaches based on inclusion of vector 
mesons into the Chiral perturbation theory framework |3i S H ■ 

No need to say that the extracted values of the meson parameters depend on the theoretical models used in this 
procedure; moreover, approximate formulas for the form factors in terms of the meson parameters inevitably introduce 
a systematic error which is extremely hard to control. 

Neglecting systematic uncertainties may lead to controversies in the determination of the resonance parameters. To 
illustrate this statement, let us turn to Tabled which presents the vector meson masses as quoted in the last three 
editions of the Particle Data Group One clearly sees a very small error and relatively sizeable 'time- variations' 

Table 1: Values of the meson masses as quoted in last three editions of PDG [lllll^. 





1998 


2000 


2002 




770.0 ±0.8 
781.94±0.12 


769.3 ±0.8 
782.57±0.12 


771.1 ±0.9 
782.57±0.12 



of the average values such that some of the results from different editions are only marginally compatible with each 
other within 3(t. Moreover, the value of the p-meson mass rup = 775.9 ± 0.5 MeV as extracted from r-decays and 
e+e^ annihilation data is well above the value of the p-meson mass as obtained by averaging all data |12|. 
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The most natural explanation of this puzzle is that the systematic errors due to reliance on theoretical models may 
be underestimated. It is clear that a reliable extraction of the resonance parameters may only be reached if as broad 
as possible a set of data and reactions is used for the analysis. Comparison with the experiment serves as a test and 
to some extent a justification of the theoretical models and approximations employed. 

The data on the e'^e~ —> 7r"'"7r~, t~ —> '!T~TT'^h'r, and e^e^ —> 7r*^7 reactions open a good possibility for extracting 
the p and uj masses and coupling constants and testing various theoretical models. 

As a first step we discuss our dispersion-theoretical framework where we start from the photon 7 and the vector 
meson V {V = p^tu, . . .) fields and identify the relevant amplitudes and couplings in a model-independent way. Then 
we introduce our model which takes into account only the most important contributions to the absorptive parts in 
the dispersion relations. In essence, (i) we keep only the contributions of the tttt and KK intermediate states and 
(ii) assume certain reduced amplitudes to be independent of the cm. energy in the resonance region. Our approach 
does not use any specific effective Lagrangian or other approximation scheme like the 1/Nc expansion. However, 
rigorous theoretical results, which are known or may become available in the future, may be easily implemented in 
our framework as improved representations for the reduced amplitudes. 

In the second step we apply our formalism to perform a simultaneous analysis of the pion electromagnetic, charged- 
current, and form factors. We obtain analytic expressions for the form factors in terms of the resonance parameters, 
paying special attention to the existing ambiguities in their definitions. We then apply our results to the Fj^, F^, 
and f-YTT data and extract masses and couplings of the vector mesons. In this paper we consider only the most recent 
data [3j Qi Us, 16] for these form factors in the range 2nijr < yfs < 1.2 GeV. We shall demonstrate that these 
data are well described by our model allowing for a reliable extraction of the vector-meson parameters. A systematic 
analysis of all available form factor data using our model is left for a future work. Then also relative normalisation 
uncertainties between various data sets will have to be considered. 

The paper is organised as follows: Section 2 gives definitions and summarises important rigorous results for Ft^ and 
F^T^. A general treatment of the vector meson-photon mixing in the framework of the dispersion approach is given 
in Section III. Our model is formulated in Section IV. In Section V numerical results are presented. Conclusions are 
given in Section VI. Appendices contain the necessary technical details. 

II. THE PION ELECTROMAGNETIC, WEAK, AND n-y FORM FACTORS 

A. The electromagnetic form factor 

The pion electromagnetic form factor is defined by 

(^+(p')|JM(0)k+(p)> =e(p'+p)^F,(g2), q = p'-p, (2.1) 

for 5^ < and by 

{TT+[p')^-[p)\J^[Qm^e[p' -p)^F^{q^), q=p'+p. (2.2) 

for q'^ > 0. Here is the electromagnetic current and e = ^/AttoZ^.- The form factor is normalised as -^^(0) — 1. 

As function of the complex variable s = q^ , the form factor i^^(s) has a cut in the complex s-plane starting at the 
two-pion threshold s — 4mJ which corresponds to two-pion intermediate states. There are also cuts related to KK 
intermediate states and multi- meson states (Stt, etc). The form factor in the timelike region (s > 0) is 

F^{s + ie)^\F^{s)\e''-^'\ (2.3) 

where 5t^{s) is the phase. For the theoretical description of the form factor in different regions of momentum transfers 
different theoretical approaches are used. 

At large spacelike momentum transfers, —q^ 00, perturbativc QCD (pQCD) gives rigorous predictions for the 
asymptotic behaviour of the form factor (17j 

F.(,^)^M^^, (2.4) 

where as is the QCD coupling parameter and /^r = 130.7 ± 0.4 MeV ^3 is the pion decay constant defined by the 
relation 



(0|J(0)7^75"(0)k+(p)) =ip^/, 



(2.5) 



3 



As the spacelike momentum transfer becomes smaller, the form factor turns out to be the result of the interplay of 
perturbative and nonperturbative QCD effects, with a strong evidence that nonperturbative QCD effects dominate 
in the region < —q^ < 10 GeV^ [isllTflj. The picture based on the concept of constituent quarks which effectively 
account for nonperturbative dynamics has proven to be efficient for the description of the form factor in this region 
(see for instance (2Qj'). 

At large timelike momentum transfers, s > 10 20 GeV^, FTr{s) can be obtained from the analytic continuation 
of the pQCD formula (|2.4|) . At small tinielike momentum transfers the situation is more complicated since there 
dynamical details of the confinement mechanism are crucial. Quarks and gluons are no longer the degrees of freedom 
of QCD leading to a simple description of the form factor. At timelike momentum transfers we are essentially in the 
region of hadronic singularities and typically one relies on methods based on hadronic degrees of freedom. In the 
region of interest to us here, — ^ 1.5 GeV^, the lightest pseudoscalar mesons are most important. There are 
many approaches to understand the behaviour of the pion form factor at these timelike momentum transfers. 

A popular approach is based on the vector meson dominance (VMD) model (jyj. In the simplest version one considers 
just the contribution of the p-meson pole, which in combination with the normalisation condition ^^(s = 0) = 1 leads 
to 



F^s)^-^. (2.6) 

TO^ — S 



This simple formula works with a reasonable accuracy both for small spacelike momentum transfers and timelike 
momentum transfers below the tttt threshold: —1 GeV^ < s < 4mJ. For s near the tttt threshold one should take into 
account effects of the virtual pions. In this region, the momenta of the intermediate pions are small and a consistent 
description of the form factor is provided by chiral perturbation theory (ChPT) 7], the effective theory for QCD at 
low energies. 

For higher s, in the region of p and lj resonances, a similar rigorous treatment of the form factor is still lacking, 
and one has to rely on model considerations. Gounaris and Sakurai (GS) _3;] obtained the expression for the p-meson 
contribution to the pion form factor which takes into account the p- meson finite width due to the virtual pions. The 
Gounaris-Sakurai form factor may be written in the form 

w^-s-^pp'^(s) 

The function -B^^(s) corresponds to the two-pion loop diagram, but one can easily add the KK loop too, see Appendix 
IbI The corresponding Feynman integral is linearly divergent, but its imaginary part is defined in a unique way. The 
real part may then be reconstructed by a doubly-subtracted dispersion representation. The Gounaris-Sakurai formula 
corresponds to the following prescription of fixing the subtraction constants 

Re B';^{s)l=^. = 0, ^Re B';^{s)l^^. - 0. (2.8) 
The phase of the GS form factor 

agrees well with the experimental data in the region 4m^ < s < 0.9 GeV^. 

The Gounaris-Sakurai form factor (|2.7|l satisfies the relation -^^^(O) = 1. But it turns out that it does not have 
enough flexibility to give at the same time a good description at the peak of the p resonance, see Appendix IbI 

Near the p-meson peak the p-meson contribution to the pion form factor can be expressed in terms of the 7 ^ p — > tttt 
matrix element as follows: 



^-(^) - T nos V (2-10) 



2.9/3— > TTTT //9 

{TT+{p')n-{p)\T\p{q,e)) = -^gp^^^ £p • [v'-pY. (0| J^(0)|p''(<z, e)) = efprupe^, (2.11 



Here ^ptttt and fp are defined according to 



where is the p-meson polarization and q is the 4-momentum vector. Now |^V(s)| from H2.10|l describes well the 
data for s ~ rn^. But extrapolating the Eq. (|2.1U|I to s = violates the normalisation condition F^(0) = 1 that is 
unacceptable. 
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Note that the definitions of gp^-mr and fp in p. 11(1 are not really appropriate since they are based on p-meson 
states. But the p meson is a resonance and has no asymptotic states. Our precise definitions of gp^T^-K and fp will be 
given below in Section IV. 

Thus, neither 1)2. 7|) nor 12.10|1 are suitable for the analysis of the form factor data for al\ s — Q~ 1.5 GeV^. There 
were many attempts to modify the vector meson dominance or to use related approaches in order to bring the results 
on the pion form factor in agreement with the data (see |^ IH 1^ and papers quoted therein) . The pion form factor 
in the region s = 1.5 GeV^ is one of the main sources for obtaining masses and coupling constants of vector 
mesons. However, with different assumptions on the form of the vector-resonance contribution to the pion form factor 
one obtains different values of masses and couplings. Therefore a consistent description of the pion form factor in 
this region in terms of the low-lying mesons (tt, p, uj) is crucial for extracting reliable values of these parameters. 
Interesting results have been obtained by the authors of who noticed that an effective momentum-dependent p7 
coupling appears in the framework of the effective Lagrangian approach. This momentum-dependent p7 coupling also 
improves the description of the pion form factor at timelike momentum transfers in the region < < 1 GeV^ . 

Clearly, to achieve a realistic description of the form factor, one has to account for vector meson finite- width and 
mixing effects. This may be done in a consistent way within a dispersion approach and will be the subject of Section 
III. 

B. The weak form factor 

The 7r~ tt" weak transition form factors parametrise the charged-current transition amplitude as follows 

(7r+(p)7r«(p')|S(0)7Md(0)|0) = V2F+{q^){p' - p)^ + V2F- {q^)q^. (2.12) 

In the limit of the exact isospin symmetry we have = and F^ — i^^|isovoctor. In practice we expect that the 
electromagnetic form factor should be close to for < < 1 GeV^ except for the region of the uj resonance. 
The form factor F^^ contains an important isospin-violating contribution of the ui resonance, whereas there is no 
contribution analogous to lo in 

C. The 717 form factor 

We shall be interested in the process e+e^ ^ 7* ^ 7r°7 where one of the photons is real and the other is virtual. 
The form factor F-^t^ relevant for this process is defined [I^, according to 

(^°(p)7(g',e)|Jp(0)|0) = e'e^f,,,e*''q^q'PF^,{q^). (2.13) 

In terms of this form factor the e+e~ — s- t:'^^ cross section reads 



2 o o / m 



,2 \ 3 



ae+e-^.o^(g2) = -TT^aL. - \F^.{q^)\'' . (2.14) 

In the chiral limit, m"^ — 0, the value of the form factor for = is fixed by the Adler-Bell-Jackiw anomaly |2^ 

i^7.(0)U2^o = TTlhrT- (2.15) 

In reality the pion is not massless, but still the anomaly provides a very good description of the observed tt'^ — > 77 
decay rate. 

We shall therefore use the value H2.15|l also for the physical pion in order to fix the form factor F^t^{Q) in our 
analysis. 

In the region of large spacelike momentum transfer g, the form factor can be calculated from pQCD with the result 



F,Aq^) - (2.16) 
Brodsky and Lepage proposed a simple formula 

^7^(9') -T-^S^, (2.17) 



47r2/^ - q 



5 



which interpolates between = and — > — cx) and works well for all < 0. This formula may be written as 




(2.18) 



with Mrcs — Ztt/tt — 880 MeV, not too far from the masses of p and u which give the dominant resonance contribution 
to the form factor. 

To describe the form factor in the region Q < q^ < 1.5 GeV^, we should again use the meson degrees of freedom. 
For a realistic description of the form factors we must take into account finite-width and meson mixing effects. 



III. MIXING OF VECTOR MESONS: PROPAGATOR MATRIX AND VERTEX FUNCTIONS 

In this section we present model-independent considerations on the mixing of the photon with vector mesons. 

Since vector mesons are unstable particles, one of the possibilities is to start with hypothetical stable states, which 
then get a width by inclusion of some interactions. This is an inherently perturbative picture which emerges for 
instance when the l/N^ expansion is used. 

We shall avoid such a perturbative approach and instead start with properly defined renormalised field operators 
with the quantum numbers of the vector mesons we are interested in. Clearly, such field operators can be defined 
in the framework of QCD. As the second step, we shall analyse the propagator matrix describing the mixing of 
these vector-meson fields with the photon field. Then, we define certain transition amplitudes (or vertex functions) 
which are one-particle irreducible in the s channel, and establish the connection between these vertex functions and 
the experimentally measured form factors. These considerations are fully general and do not include any model 
assumptions. As the next step our model is formulated making certain assumptions for these vertex functions. This 
procedure is similar to the one used in in the discussion of the "f-J /tp mixing. 

Let us consider the photon field A^{x) and a set of hermitian neutral vector-meson fields vll^\x), (j = 2, ...,n). 

For convenience of notation we set vll'^\x) = A^{x). 

The fields vli'\x) (j = 1, ...,n) have the same quantum numbers and therefore will have anx n propagator matrix 
describing their mixing 

^^^"^(9) = 7/ d'a;e'^-(0|r*{Fj^)(a:)T/W(0)}|0). (3.1) 

Here T* is the covariant version of the T product, see for instance Using a covariant gauge for the photon, we 
can separate A^f^'^'' into transverse and longitudinal parts as follows 

Aif;'^) = (-9,. + 4^) A^^-'\q') - 4f- ^f^i^)- (3.2) 
^ \ q- + q- + ie 

The matrices 

At,l(<z^) = (A^i^g^)) (3.3) 

are analytic in the complex q^-plane with cuts on the positive real axis. We shall always work to leading order in 
the electromagnetic interaction. Then the leftmost cut starts at q^ — Arn^, the two-pion threshold. Similarly, the 
transverse and the longitudinal structures can be isolated in the inverse propagator matrix 

(A-(.))i'=^ = (-.9.. + f^) {A-Al')f-'' ~ ^ {Al\l')f-'' . (3.4) 

The propagator matrix satisfies several general relations. 
• Translation invariance of the vacuum gives 

A(i'=)(9)-A(^^)(-g). (3.5) 



• CPT invariance gives 

A(/^'=)(g) = A(f-'^)(-g). 



(3.6) 
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• T-invariance of strong and electromagnetic interactions gives 

A(f;'=)(<?°,<7) = A(^''=)'-(-g«,<7). (3.7) 

From (|3.5|l we find that the matrices At,l must be symmetric 

{^TAe)f = AT,L(g'), (3.8) 

whereas ()3.6|l and (|3.7() are satisfied automatically and give no restrictions. 

Using next the hermiticity of the fields V^i''^ (x) we get the unitarity relation for the propagator matrix 

A(/''^-)(9)-(A(^^^)(q))* = -^| d'xe^''^{0\{v^^\x)V}'\0) + v}'H0)V,9H^)}\0) 

X 

+ {2iTf6^'\q+pxmv!;''\nx{px)){X{px)\vli=\m)] (3.9) 

Here we have inserted a complete set of asymptotic (in strong interactions) states \X{px))^ where px is the four- 
momentum. Note that the states \X{px)) contain pions and kaons, but no p or w mesons since the latter are unstable 
and thus have no asymptotic states. 

Let us now define for all states \X{px)) the reduced, or amputated, matrix elements (X(px)||V';i"'''||0) by taking 
out of {X(px)\V^^'^\Q) all s channel ^-propagator terms: 

{X{px)\\V^/^m = {X{px)\V^^>\Q) (A'i(px))l^' . (3.10) 

Here and in the following we use the summation convention. The inverse of (13.10(1 reads 

{X{px)\Vli^M - (X(px)||^«1|0)A(;f)(p;,). (3.11) 

The reduced matrix elements, or vertex functions, (X(px)||(Vji"'''||0) are one-T^ irreducible in the s-channel. 
It is convenient to define the transverse and the longitudinal components of the vertex functions 

(^bx)||T4?||0) = (X(px)||W^>||0)(5..-^^^^^ 

{x{px)\\vi^'^m = {x{px)\\v^'>m^^^. 

Now insert (|3.11() into the unitarity relation (|3.9|l . Considering first (|3.9() for > 0, we obtain 

A^^^Hq) - Uill'Hq))* - -^^(2^)^<5(4)(g ~p^) (aJ;^))* {Xipx)\\V'-^'^^'m^X{px)\\V^''>'mAi''f\ 



(3.12) 



X 



(3.13) 

Multiplying H3.13|l by (A^^)^ from the left and by A^^ from the right we come to the unitarity relation for the inverse 
propagator 1(3. 4|l as 



where the discontinuity matrices Dt.l{q^) are given by 

X 

D^t'\q') = -i$^(2^)^<5(4)((7-p^)(X(p^)||F«||0)*(X(p;,)||Ff)||0). (3.16) 



2 

X 
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and 



The discontinuity matrices satisfy the relations 

Driq^) = 0, for q^ < Ami, 

Driq^) > 0, for q^ > Ami, (3-17) 

DUq') - {DUq')f = {DL{q'))\ 

DL{q^) = 0, for q2 < 4^2^ 

Ddq^) < 0, for q^ > Ami. (3-18) 

Considering in 1)3.9(1 the case (7° < 0, inserting 1(3.11(1 and using 1(3.5(1 we find exactly the same relations 1(3.14(1 - 1(3.18(1 . 

In our applications the longitudinal part A^^q^) plays no role, so we concentrate on the transverse part Axiq^) in 
the following. 

The analytic properties of A^^[q'^) allow us to write a dispersion relation for it, which we assume to be convergent 
with two subtractions: 

A^\q^)^-M' + Kq^ + iqY^ f ds ^^^f> (3.19) 

■K J s'^ [s ~ q'^ — le) 

Here the subtraction terms Af^ and K have to be constant real symmetric matrices 

= (M^)*, K = K'^ = K*. (3.20) 

This is as far as we can come with a general analysis of the propagator matrix. 

In the next section we shall analyse the amplitude {X{px)\Afi\Q) for the electromagnetic field and the state \X(px)) 
being the tttt state, which gives the pion form factor. The equation ((3.11() with j = 1 represents this amplitude in 
terms of the vertex functions 1 10) and the propagator matrix Ap,y for which we have the dispersion 

representation. 

The merit of this representation is that different types of singularities are isolated in different quantities: the 
propagator matrix contains the resonance poles which lead to " fast" variations of the form factors in the resonance 
region; the reduced amplitudes are free from these singularities and therefore represent slowly varying functions in 
the resonance region. To go further with the form factors we need some dynamical inputs for the vertex functions 
{X{px)\\V^^^''\\0) and for the matrices and K, see Section IV. 

Before going to the details of the model, note that we are free to change the basis for the fields. Defining new fields 

V;(^)(x) = QfeV;W(x), (3.21) 
with C = {Cjk) a real non-singular n x n matrix,^ we get the propagator matrix of the new fields 

A^,{q) = CA^,{q)C^. (3.22) 

This leads to 

At,l(9') = {C-')"^ A^^^iq')C-\ (3.23) 

(^(^'^)I|V;(^'||0) = {XipxW^'^mC^l. (3.24) 



^ In field theory we have also a freedom to make more complicated redefinitions of the fields, for instance 

Vij.{x) = {1 + cn}Vfx{x), c = const. 

Such transformations will change the behavior of the propagators and the behavior of the vertex functions. In the present article 
we will not explore further the possibility of such field redefinitions. 
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The freedom of the field redefinition (|3.21() can and will be used to impose certain constraints on the matrices and 
K in 13.19|l . If if is a positive-definite matrix - as it should be from the positivity of the metric for physical states 
in the Hilbert space - we can, for instance, diagonalise K and simultaneously by a transformation (|3.2HI . The 
procedure to achieve this is completely analogous to the introduction of normal coordinates in the problem of small 
oscillations around a stable minimum of the potential in mechanics (see for instance 33]). 

We should, however, be careful with redefinitions of the photon field vll^^ — A^. A redefined photon field containing 

components proportional to the strong interacion vector fields V/i"''' with j > 1 will induce a direct quark-lepton 
coupling. We think this is unacceptable. The conditions which allow to avoid this and to guarantee the massless 
photon and the correct charge normalisation are summarised in Appendix^ and lead to 

Ml = 0, j = l,...n, (3.25) 
Kii = 1. (3.26) 

IV. THE J- p-Lo SYSTEM 
A. The model 

We now apply the general considerations of the previous section to the system containing the photon field and the 
vector meson p and uj fields, 

VW(a:)^A,ix), V(^\x)^p,{x), V^^\x)^Lo,{x). (4.1) 

We suppose the field to be purely isovector, and lo^ to be purely isoscalar. The electromagnetic coupling and isospin 
breaking from different up and down quark masses in QCD will introduce non-diagonal terms in the propagator matrix. 
In the following we will frequently use the indices 7, p, and u instead of 1, 2, 3. 

The 3x3 matrix of (|3.19|) . (|3.2U|) for our system has to satisfy H3.25|l . By a linear transformation H3.21|l . but 
involving only the p and lo fields we can make diagonal. 

At this stage we define the p and oj mass squared parameters and as zero points of the real parts of the 
diagonal terms of the inverse propagator matrix, that is by the relations 

Re ( 1) [ml) ^ 0, Re {^^') (m^) = 0. (4.2) 

Then we choose the normalisation of the fields and uj^ such that the matrix has the form 



/ \ 

. (4.3) 

yO ml) 



To calculate the dispersive part of the inverse propagator, we must restrict the set of the intermediate states \X{px)) 

to be included in the unitarity relation H3.13|l . and parametrise the reduced amplitudes of the fields V^^'' between 
these states and the vacuum. 

Assumption 1: as the intermediate states \X{px)) in the dispersion relation (|3.19(l with Dt given by 13.15(1 we 
shall consider only tt+tt", 37r, K^K", and K^K^ states. 

For the tt+tt" and KK states we have 

(7r+(fci)^-(fc2)||V^J|0) = gi^l{k,-k,),, 
(X+(fci)/f-(fc2)||y^J|0) - 5^)^(fci - (4-4) 

where g'^l and g^^j^ are in general (slowly varying) functions of (fci -I- A:2)^- 

Assumption 2: in the region of interest we neglect the dependence of g'^Jl and g^xK '-''^ (^1 + ^2)^ and assume all 
g^') to be real constants 

„(i) _ „ „(2) _ i„ „(3) _ 1 u r.\ 
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and similarly for the KK intermediate states 

(1) (1) A (2) (2) 1 (3) (3) 1 / . „x 

9k+K-=^' 3x0ifo=0, 9k+K- = ^9kok'> = 29P^KK, gK+K- = Sk^^ko ^ l^gi^^KK- (4.6) 

Here g^^} — g^jt\r^_ — e and g\)ai?o = as required by the charge normalisation of the 7r+, and if*^, see Appendix 

The decays p — > KK and lo — > KK are forbidden kinematically at the p and lo peaks. This makes a direct 
determination of the corresponding coupling constants Qp^kk and g^j^KK difficult. Therefore we use as additional 
theoretical input the relations following from the approximate SU(3) flavour symmetry of strong interactions and 
ideal mixing of the vector mesons 

gp~*KK = gu-^KK = -^gp^TTTr- (4.7) 

In this paper we do not analyse the Stt decays in detail. In the unitarity relation (|3.13() , (|3.14() the 37r intermediate 
states produce the width of the lu, T^, which is one of the fitting parameters. 

We now have to specify the matrix K in H3.19|l . The explicit form of the matrix K is discussed in Appendix^ and 
here we present the final form for the inverse propagator matrix in our model: 



/ s e^s + B^p{s) e^s + B^^{s) 



e^s + -S7p(s) -ml + s + Bpp{s) sbp^ + Bp^{s) \. (4.8) 

V e^^s + B^i^(s) shp^ + Bp^{s) -ml + s + B^^{s) 



The functions Bij are constructed by the doubly-subtracted dispersion integrals (|3.19|) corresponding to the pion and 
kaon contributions and include also the relevant subtraction terms defined such that 

By (s = 0) = 0, ReBpp{ml)^ReB^p{ml)^ReBpUml)^0, ReB^^iml) ^ ReB^^{ml) ^ 0. (4.9) 

For the functions Bij defined according to the conditions (|4.9I) . the dimensionful constants fp and correspond to 
our precise definitions of the leptonic decay constants of the vector mesons. The detailed formulas for Bij are given 
in Appendix IbI 

The intermediate tt^ states do not contribute to the form factors to first order in the e.m. coupling. Nevertheless, 

we need the reduced ttj amplitudes for the description of the ttj transition form factor. The reduced amplitudes 
have the form 

(7r(fci)7(fc2,e)||T^«(0)||0) = ee^^^.e-fc^^f , (4.10) 

where g^} in general depend on (fci + fc2)^. We assume the g^} to be constant as we did for the tttt and KK couplings. 
Then g^^J is determined from the anomaly (see appendix C) 

.«=e^,.(0) = e^^. (4.11) 
The two additional dimensionless parameters Qp^-y-K = Q^-rr and g^^j^^- 

— ttIqj Q'jtt are assumed to be real. 

Let us summarise the parameters of our model. These are: 

• the Breit-Wigner masses of the vector mesons, that is rrip and m^j, 

• the decay constants fp, f^^, 

• the mixing parameter &p^, 

• the couplings of the vector mesons p and u> to two pions gp^-KTr, guj^T^r, 

• the width of the uj meson T^, 

• and the 717 couplings of the vector mesons gp — ^7^7 and g^^ — ^777- 

Values for these parameters can be found by the fit to the available form factor data. However, it turned out that the 
parameters fuj and can not be well determined in this way from the reactions under discussion. The right place 
to extract these parameters from the experimental data is the reaction e^e" — > 7* — > 37r, where the contribution of 
u! dominates. We leave a study of this reaction for a separate paper. Here we fix /^^ and to the PDG values, and 
leave only the remaining parameters from the list above free in the fits. 

Clearly, the inclusion of higher resonances with the same quantum numbers in the mixing scheme is straightforward. 
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B. The form factors 



The calculation of the form factors is now straightforward: we must reconstruct the propagator matrix from its 
inverse, and then calculate the amplitude (X|A^|0) from (|3.11|) for the relevant final states. Finally we have to take 
into account that the amplitude of the e.m. current {X{q)\Jp\0) is related to the amplitude of the electromagnetic 
field as 

{Xiq)\Jp\0) = {-gp,q^ + q^q,) {X {q)\A''\0) . (4.12) 

For the pion form factor the final state is the tttt state, and for the ir^ form factor it is the ttj state. Since we work to 
first order in the e.m. coupling, the set of the intermediate states is the same in both cases. In the model described 
above it includes the tttt, Stt, and KK intermediate states, and we use (|4.4|l and (|4.10() for the vertex functions. To 
first order in the electromagnetic coupling the expressions for the pion elastic and the 7r7 form factors obtained by 
the procedure described above may be written in a simple form (see Appendix 

F^is) - l-sG^^,,A(s)G^^y, (4.13) 
F^^is) = F^^{0)-.sGl._,,^A{s)G^^v, (4.14) 

with FjTr(0) given by (|4.11|l . The propagator matrix A here is the inverse of the vector meson block of the matrix 

A-i (ESI) 

^-1 ^ 1 + s + Bpp{s) sbp^ + Bp^{s) I 

V sbp^ + Bp^{s) -ml + s + ' ' ' 

and 




9p- 




Gv^.. = I f^"^"" 1 , Gy^^, = ^ j , G,^v = I 7; , ) . (4.16) 

In the case of the charged-current form factor F+ (I2.12() describing the tt" — ^ 7r° transition the w contribution is 
absent so we find from (|C16p of Appendix lO 

m2 - s - Bpp{s) 

Note that in the expressions for the pion electromagnetic and 717 form factors H4.13|l and (|4.14|) . the parameters rup 
and fp are those of the p*^ meson, whereas in the formula for the charged-current form factor H4.17|l the parameters 
refer to the p~ meson. 

The expression ()4.17|l can be written in the "usual" vector meson dominance form H2.10|) 

n^l-s-Bppis)^ (4-1^) 

in terms of an effective s-dependent 7/3 coupling 

„ s 2(m^ — s) 
ffis) = fp^ + ^ '- (4.19) 

9P — ^TTTT^p 

In this way we can make contact with the results of Q where an s-dependent 7p coupling is defined in an effective 
Lagrangian approach. 

We use the expressions (|4.14f) and (|4.17|) for the numerical analysis of the data for the pion electromagnetic, charge 
current, and ttj transition form factors in the next Section. 
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Table 2: Parameters of the resonances as found by the fit to the form factors in the region ^/s < 0.9 GeV within the "f — p — uj 
mixing scheme (Fit I) . Higher resonances are not included. The PDG value for is used. Fit to and F^^: x^/DOF = 75/74. 
Fit to the charged-current form factor F^: /DOF = 11/23. The pui mixing parameter has the value hp,^ = (3.5 ± 0.6)10"^. 



Res. 


mv, MeV 


fv, MeV 






P~ 


775.5 ± 0.4 


152.5±0.33 




11.52± 0.04 


P° 


773.6 ±0.5 


154.1±0.67 


0.60± 0.06 


11.43± 0.04 




782.42 ± 0.04 


45.3 ±0.9 


1.79± 0.09 


-0.27± 0.13 



Table 3: Parameters of the resonances as found by the fit to the form factors in the region -^s < 0.9 GeV (Fit II). The 7 — p — lj 
mixing scheme with addition of the p' = p(1450) is employed. The PDG values for fuj and mpi are used. Fit to Ftt and F-,Tr: 
X^/DOF = 68/72. Fit to F+: x^/DOF = 11/21. The parameters of p' cannot be determined by this fit. The pco mixing 
parameter has the value bpuj = (3.7 ± 0.6)10"^. 



Res. 


mv, MeV 


fv, MeV 


gV^TT-y 


gV^TTTT 


P~ 


775.3 ± 0.8 


152.4±0.4 




11.50± 0.05 


p' 


773.8 ± 0.6 


155.3±3.2 


0.61± 0.06 


11.53± 0.10 


LU 


782.43 ± 0.05 


45.3 ±0.9 


1.76± 0.09 


-0.31± 0.10 


P' 


1465 ± 25 









Table 4: Fit to the data for < 1.2 GeV (Fit III), where p, u, p' = p(1450), and p" = p(1700) are taken into account. Fit 
to F^ and F^^: x^ jDOF = 72/89. Fit to F+: x^ jDOF = 13/27. The extracted p and lo couplings and masses are very stable 
with respect to inclusion/exclusion of p". The couplings of p' and p" correlate very strongly and cannot be reliably determined 
by this fit. The piu mixing parameter has the value hp^ = (3.5 ± 0.5)10"^. 



Res. 


mv, MeV 


fv, MeV 


gv^TT-i 


gV-^TTTT 


P~ 


775.3 ±0.5 


151.5±1.5 




11.50± 0.05 


P° 


773.7 ±0.4 


155.4±1.7 


0.65±0.05 


11.51± 0.07 


u; 


782.43 ± 0.05 


45.3 ±0.9 


1.73±0.08 


-0.35 ±0.10 


p 


1465 ± 25 








p" 


1700 ± 20 









V. NUMERICAL RESULTS 

We analyse the recent data on the pion electromagnetic form factor H2.2|l . the charged-current form factor 
(|2.12|) . and F^^, form factor H2.13|l in the region = 1.2 GeV using the formulae H4.18|l . H4.14|l . and H4.17|l . 
respectively. We take into account that the pion electromagnetic form factor and the ir^ transition form factor 
-FIyt contain contributions of the neutral p° and w resonances, whereas the charged-current form factor F^ contains 
the contribution of the p~ meson. Since we consider the isospin-violating — lo mixing effects, we do not assume 
the parameters of the charged and the neutral p mesons to be equal to each other. We therefore fit the data for the 
F.^. and F.^^ form factors and extract in this way the w and (P parameters. We use the recent SND data for the 
form factor F^t^ and the recent update Ti'l of the CMD-2 data '15*1 for F^r. We also include the available data on the 
phase of the electromagnetic form factor 30] . 

We separately fit the form factor F+ and extract the p^ parameters from the CLEO data ■ 
We perform different fitting procedures explained below. The fitted values of our parameters are given in Tables 
EEl andH 

The first fitting procedure (Fit I) includes the p^ , and uj contributions, and neglects effects of higher resonances. 
We work in the region < 0.9 GeV where this approximation is checked to be self-consistent. The x^ /DOF for 
this fit is 75/74 for F^- and F-y^. It is 11/23 for F+. The resulting values of the parameters are given in Table |21 

In the second step (Fit II), we study the stabihty of the extracted p and w parameters with respect to the inclusion 
of higher resonances. This might serve as a probe of the systematic errors. We still stay in the region ^/s < 0.9 GeV, 
but include in addition to p and w also the p' = p(1450) resonance, assuming its contribution to have the form (I2.10f) 
with the corresponding mass and coupling constant of p' . We fix the mass of this resonance according to PDG 12] and 
neglect possible p — p' mixing effects. The parameters obtained by this procedure are given in Table |3| The quality 
of the fit to Fn- and F^^ definitely improves (68/72), whereas for F+ no improvement is seen. It is worth noting that 
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Table 5: The pion form factor at small momentum transfers. Coefficients of the expansion 15.111 from fit III are given. The 
errors given in our results are only statistical emerging from errors in masses and couplings. 







cl, GeV"* 


Our result 
ChPT to order 0(/) [35] 
Results from [9] 
Naive VMD Eq. JTSJ 


11.41±0.05 
11.22±0.41 
11.04±0.3 
10.16 


3.83±0.02 
3.85±0.6 
3.79±0.04 
2.8 



the extracted parameters of the p and lo turn out to be very stable with respect to the inclusion/exclusion of p'. 

In the third step, we extend our analysis to the region ^/s < 1.2 GeV (Fit III). We fit the form factors taking 
into account the three resonances p, uj, and p' {x^/DOF = 72/89 for and F^^; x^/DOF = 13/27 for F+). The 
results of this fit are given in Tabled We then also include in addition to them the p" ~ /3(1700). The coupling 
constants of p' and p" turn out to be strongly correlated with each other and therefore cannot be extracted from 
the data under consideration. For a reliable extraction of these parameters one should go to higher values of s. 
Important for our analysis is that the masses and couplings of the p and lo mesons are remarkably stable with respect 
to inclusion/exclusion of p" and very well compatible with the numbers obtained in Fits I and II. The form factors 
calculated with the parameters from Table 01 are shown in Figs^l^l and 13 as solid lines. 

The very satisfactory description of the data speaks in favour of the reliability of our assumptions of the dominance 
of the TTTT, KK, and 37r intermediate states and on the negligible s-dependence of the vertex functions. 

Note that the masses of the charged and the neutral p mesons are different as obtained by our fits exposing an 
isospin violation in p mesons which is extensively discussed in the literature (see [33 | and refs therein). 

However we would like to point out that assuming the parameters of the charged and the neutral p mesons to be 
equal to each other also leads to a very good description of the data with /DOF below 1 for all fitting procedures 
l-III (Fit I: x^/DOF = 94/100, Fit II: x^/DOF = 83/98, Fit III: x^jBOF = 86/121). Therefore, strictly speaking, 
the data analysed by us here does not require the masses and couplings of the charged and neutral p mesons to be 
different. Still, we do not think it reasonable to take into account isospin-violating p ~ uj mixing effects and assume 
the absence of such effects in the charged and neutral p mesons. We therefore do not discuss in detail the fitting 
procedures in which the parameters of the charged and neutral mesons are assumed to be equal. 

Fig. 21 shows the pion form factor at small spacelike momentum transfers, which was not included in the fit. Still 
one can see a good agreement with the data even down to s = — 2 GeV^. The improvement of the description of F^^ 
compared to the naive VMD Ansatz H2.6|l shown as dashed line is obvious. Similar results were obtained in |5J. 

Next we study the low energy expansion of the pion electromagnetic form factor near s = 

i^.(5)-l + i(r2)^ + 4^2 + 0(^3). (5.1) 

Here {r'^)y is the squared charge radius of the pion. This quantity and Cy are of great interest in the framework of 
ChPT for fixing certain parameters, see for instance 9, 35]. In Table Owe compare our results with those of 9, 3^ 
and the naive VMD. We find full consistency with ChPT within the errors which for our results are only statistical 
ones. Of course a study of systematic errors should also be done but this is beyond the scope of this paper. 



Finally, let us discuss the p and w masses and the p width. The value of the p-meson Brcit-Wigner mass defined 
in (|4.2|l nip = 774.3 ± 0.4 MeV agrees with the value obtained recently from the weak pion form factor jl^], and is 
sizeably higher than the value rup = 771.1 ± 0.9 MeV quoted by PDG. 

The Breit-Wigner width of the p meson is defined according to the relation |lO|| 

l/Tp ^rup \s=ml, (5.2) 

where (5^(s) is the phase of the pion form factor F^{s). Numerically, we obtain from fit HI for the charged p 

F^^ = 149.85 ±0.4 MeV. 

Next we turn to the pole masses and widths of the p^ , p'~' and uj. For the charged p-meson the location of the pole 
in the second Riemann sheet of the s-plane is found by solving the equation 



ml - s - Bpp{s) = 0. 



(5.3) 
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Table 6: Pole masses and widths of p and u. 











A/P°'=, MeV 


756.7±0.4 


755.8±0.4 


782.44±0.05 


^P°'^ MeV 


144.7±0.4 


143.8±0.4 


8.38±0.05 



The corresponding solution, Spoie, can be used in two different ways to define pole masses and pole widths. Either we 
choose to set 

Spoie = M'l - zr>i<=M'p, (5.4) 

or 

spoic = (Mp - iTl°'-l2f. (5.5) 

The values obtained with the definition H5.5|l are given in Table El 

The pole masses of the p° and the w meson are affected by the poj mixing effects and are obtained from the equation, 
see ^ 

{ml-s-Bpp{s)] {ml-s-B^^{s)]-{sbp^ + Bp^{s)Y ^G. (5.6) 
The corresponding values of the pole masses of the and oj and widths are also given in Table 



VI. CONCLUSIONS 

We have discussed a general approach to the description of vector mesons and their mixing using dispersion relations. 
This approach allows us to represent various observables for the vector mesons as products of the vector-meson 
propagator matrix and the reduced amplitudes (vertex functions). The unitarity relation gives the anti-hermitian 
part of the propagator matrix in terms of the relevant reduced amplitudes. 

The merit of this approach lies in the possibility to separate different types of singularities in different quantities: 
the propagator matrix of the vector- meson fields contains the resonance poles which lead to "fast" variations of 
the form factors in the resonance region, whereas the reduced amplitudes are free from these singularities and are 
therefore slowly varying functions in the resonance region. The description is fully general at this stage and contains 
no approximations. To go further we need some dynamical inputs for the reduced amplitudes. 

We then formulate our model for the form factors based on the following assumptions: 

(i) we take into consideration the resonances p, lo, and in a rough way also the /9(1450) and neglect higher vector 
mesons; 

(ii) we take into account the tttt, KK and effectively also Stt intermediate states, and neglect contributions of 
multi-meson states in the unitarity relations; 

(iii) we assume the scalar coupling factors in the reduced amplitudes to be constant in the region of the momentum 
transfer = - 1.2 GcV. 

On the basis of these assumptions we perform a combined analysis of the recent data for several reactions: e^e^ — s- 
7r+7r~, T~ TT^n^Vr, and e+e~ in the region of ^/s = — 1.2 GeV. All the analysed data is well described in 

our approach, allowing for an extraction of the resonance parameters, such as the Breit-Wigner masses and effective 
coupling constants. Our main numerical results are given in Tables EBl 

These results are obtained by the fitting procedures, which allow for different masses and couplings of the charged 
and the neutral p mesons. Therefore, we have fitted separately the charged current form factor and the neutral 
current form factors F^, and F^^r • Still we would like to point out that assuming the equality of the parameters of the 
charged and the neutral p mesons and fitting all three form factors F^ , F^, and F^t^ simultaneously also leads to a 
good description of the data with 'x^ j DOF below unity. However since the p — oj mixing and the isospin violation in 
p mesons are both effects of the same nature, we do not think it reasonable to include only the first of these effects. 
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600 800 1000 1200 



is, MeV 




600 800 1000 . 1200 

is, MeV 



Fig . 1: The pion electromagnetic form factor F^{s): (a) Modulus squared: data from li4 | ( emptvl used in the fit. Data from 
|29| (full) which was not used for the fit is shown for comparison, (b) Phase: data from |3Q]. Solid lines - the full form factor 
as obtained by Fit III, dashed lines - the p-contribution in Eq. H4.18|l . 
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-3 -2 -1 

s, GeV 

Fig. 4: The pion electromagnetic form factor: \F^{s)\^ for s < 0. Data from ,32] . Solid line: the full form factor from fit III, 
dashed line: F-^ = 1/(1 — s/m^). 



We therefore do not discuss in detail the hypothesis that the masses and couplings of charged and neutral p mesons 
are equal. 

The small errors of our results are statistical errors only indicating a good description of the data by our form 
factor formulae. Still our form factors are based on certain model assumptions, therefore a systematic error should 
be added. 

One way to estimate the systematic error is to vary the fit range, to choose different parametrizations for the 
reduced amplitudes, and to study the corresponding variations of the fitted quantities. Partially, we did this by 
choosing various ranges of s and including higher resonances in our fits. We observed a good stability of the extracted 
resonance parameters and therefore do not expect the systematic error to be large. A more detailed study of the 
systematic errors is left for future work. 

Further topics to be studied in the future are the inclusion of the cj) meson in the mixing scheme and the implications 
of our form factor formulae for the theoretical values of the .g — 2 factor of the muon and the fine structure constant 
at the Z mass, as{mz). 

To summarise, we have obtained a very good description of the form factors in a model based on the unitarity 
and dispersion relations. Chiral perturbation theory constraints were checked to be respected. We have introduced 
reduced amlitudes or vertex functions which describe the coupling of the photon 7 and of the vector meson p and 
u fields to two pseudoscalar mesons and to n^. These vertex functions contain invariant coupling functions which 
in principle depend on the momentum transfer ^/s. Perhaps the most interesting result of our study is that these 
invariant coupling functions are really coupling constants for < -^/s < 1.0 GeV. Thus these couplings are "frozen" 
below the GeV scale. This may perhaps be related to the expectation mentioned frequently in the literature |,3fi,] that 
in QCD the coupling parameter as{p'^) is "frozen" for p^ ^ 1 GeV^. 
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APPENDIX A: GAUGE INVARIANCE AND PROPERTIES OF THE MATRICES AND K 



Here we study the conditions to be imposed on the propagator matrix H3.1|l arising from the requirements to have 
a massless photon, the correct charge normahsation, and no strong-interaction long-range force. This leads to the 
following constraints at 0: 



1 



= _ (1 + 0(^2))^ 
Ai^^'Hq') = 0(1), for(z,j)^(l,l). 



These conditions require the matrices M and K to satisfy 

r2 



Mt^ = 0, forj = l,...,n, 
Kii = 1. 



(Al) 
(A2) 

(A3) 
(A4) 



To show this, we note first that by a transformation of the type ()3.21|) . but involving only v/^''^ with j > 1, we can 
always achieve 

Then the matrix takes the form 



0, forj>2. 



(A5) 





V 



Ml, 



0\ 



(A6) 



where 



For n > 2 we define 



Mh ■ 


■ ML 


Ml, . 


■ ML 


Mis ■ 


■■ Mi^ 


Ml, . 


■■ ML 



(A7) 



(AS) 



and for n = 2 we set AP = 1. 

In the following we assume that Af^ is a positive-definite matrix. This means, for instance, that all fileds 
j — 2, . . . ,n must be independent. 

From H3.19|l and (|A6|I we find now for ^ 



det Ay^(g2^ 



= {-iy'det]VP + OiqL 

= (-1)" [Mi^detM^' - (M^^fdetM^" 



(A9) 



(14). 2, _ det {-hp ) + 0{q^) 



detAj;i(g2) 
-det M^' -f 0(^2) 
detM2-hO(g2) ' 



(AlO) 



(1,2), 2, _ (^l)"M2,detM2"+0(g2) 



detAj;i(g2) 
M2^detM2" -H0(g2) 
detM2-HO(g2) ' 



(All) 
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We see now from (|A10ll that in order to fulfil IjAip we must have det = 0. Then HA2|) for = (1, 2) and HA11|) 
require = 0- Combining this result with l|A9|l gives M^^ = 0. Recall that we assumed > which implies 
also > 0. Thus, we have already demonstrated (|A3|) . Inserting now (|A3p into H3.19|l gives 



With (fnol) this leads to 



detAj,\q^) = (-l)"-ig2^iidetM2' + 0((g2)2). (ai2) 
(-l)"-idet M2' +0(92) 



To fulfil (|A1|) we must, therefore, have Kn — 1 which proves l|A4|l . 



(-l)"-ig2Xndet M^' + O ((g2)2) 

^ + (A13) 



APPENDIX B: THE INVERSE PROPAGATOR MATRIX FOR THE j - p - uj SYSTEM 
Given H4.3|) , the inverse propagator (|3.19|) for the ^ ~ p ^ uj system can be written as 





ml, I iml 



^t\s)^-\ ™? I +sK + s^- J ds^ ,,^,^^^^^ . (Bl) 



Here we set = s and Dt{s') is obtained using Assumptions 1 and 2 of Section IV and Eq. H3.15|l . 
It is convenient to split the constant matrix K into two parts 

K^K'' + K^, (B2) 

to be specified later, and to define the matrix function 

00 

B{s) ^ sK^ + s^^- f ds' .. . (B3) 

TT J s''^[s' ~ s ~ le) 

We thus have 

l-{Bis)-BHs)}^DTis). (B4) 

The matrix if will be chosen such that it cancels the real part of the dispersive contribution to B{s) in Eq l|B3l) at 
some specific values of s. Then splitting the matrix K according to ljB2|l has an unambiguous physical meaning. 
We choose -fC"^ = 0. With this, the condition (|3.26ll requires if^^ = 1. 

The diagonal elements Kyy, V = p,uj, are chosen such that they cancel the real part of the dispersive contribution 
in ljB3|l at s = my for the vector meson V. That is, we require 

ReBppiml) = 0, ReB^^{ml) ^ 0. (B5) 

The normalisations for the p and u fields are fixed by requiring (|4.3|l . Given ljB5|) . this requirement leads to 

kL - = 1. (B6) 



The non-diagonal matrix elements {K'^)jv are chosen to cancel the real part of the dispersive contribution to B^v 
2 

V 



at the point s = my. That is, we require 



Now we set 



Re B^p{ml) = 0, Re B^^{jnl) = 0. (B7) 



{K%v = e^, (B8) 
mv 
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where the real parameters fv correspond to our precise definitions of the leptonic decay constants of the vector mesons 
V = p,uj. 

Finally, the matrix element {K°')pu} is chosen such that it cancels the real part of the dispersive contribution to Bp^^ 
in ljB3p at the point s = ni^. Thus we require 



Then the element 



ReBp^iml) = 0. 



(B9) 



(BIO) 



defines our p—uj mixing parameter. Clearly the value of the mixing parameter depends on the choice of the subtraction 
point in (|B9|) . but physics, of course, does not. 

Let us now discuss in more detail each of the functions i?.y(s) of (|B3|1 . 



a. B„ 



The function Im Bpp receives contributions from the tt+tt , K~^K and K'^K'^ intermediate states. The two-pion 



■•pp 

contribution to lia Bpp reads 



where 



lmBpp{s)\^^ = 
It[iB^^{s) = I{s,ml), 

J(s, TO^) 



(Bll) 



1 / 4m2\3/2 

s i ' '^-^^^ 



1927r 



(B12) 



We have to take into account also K^K~ and K^K° intermediate states which give contributions similar to I|B11|I . 
(|B12p with TT replaced by K. The coupling constant gp^KK cannot be measured directly so we use the flavour-SU(3) 
relations H4.7|l 



Then we find 

l-mBpp{s) = g 
As explained above, see (|B4|I . we require 



,2 



'^gp^KK — 9p- 



Im Btttt + ^ (Im Bx+K- + Im -B^o^o) 



Re Bppim^) = 0. 



Putting everything together we find from (|B3|I 



y p — >7r7r 



Im +2^"^ 



(B13) 



(B14) 



(B15) 



R{s,ml) - R{ml,ml) + 



Bpp{s) = ffp^^^ 
Here I{s,m?) is defined in HB12|) . and 



^ V.P. r 

1927r2 J^„^, {s' - s)s' 



Am' 



^(i+e^ + ^log(l-| 
96^^(1 -e2 + e3arctan(i) 

where V.P. means the principle value. 



^9o 



lis, ml) + 



I{s,m\) 



= Jl 



4m ^ 



4m ^ 



- 1, 



= ^ I- 



4m^ 



for s > 4rn2^ 
for < s < 4m2, 
for s < 0, 



(B16) 



(B17) 
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b. B^^ 

When considering the B^^^^, the intermediate K^K^ , K'^K^ , and Stt states should be taken into account. 

First, recaU that the couphng constant g^^^-K is much smaUer than the couphng constant guj^xK- This becomes 
clear by using the SU(3) relation g^^^KK = ^5p-+irir- Therefore, the contribution of the Stt intermediate states can be 
safely neglected in the real part of B^^t^ . 

However, the decay uj — > KK is forbidden kinematically. Due to this, the uj width emerges mainly due to the Stt 
intermediate states and is small because guj^si^ is small. Therefore, the Stt states are crucial for the calculation of the 
imaginary part of B^^^ below the KK threshold. Still, the s-dependence of Im B^^ can in practice be neglected, and 
we may therefore use for it a constant width approximation F^^ = const. Taking into account that Bi^,^ = for s = 
()4.9|l . and making use of (|B5|I . we can write with sufficient accuracy for our purposes 

Bujuj{s) = gl^KK s [^(s, mji) - R{ml, m|.)] + iT^s/m^. (B18) 



The function Im _B^p receives contributions from the tt+tt and K^K intermediate states. Taking into account 
the SU(3) relations H4.7|l between the coupling constants we find 



Im Bjp = 2 e gp^^^r ( Im Btttt + Bk+k- 1 = 2 e ^p^tttt ( Ini Stttt + ^kk ) , (B19) 



and hence from (IBl 



2 e 

Im B-^p{s) = Im Bpp{s). (B20) 



With the conditions (jB4|l and (jB7|l we get also 



B,p{s) = -^Bpp{s). (B21) 

gp^-KTT 



The function Im B^^ receives contributions from the K~^K intermediate states: 



ImB^i^ = 2eg^^KKim Bji+K- . (B22) 
Requiring (|B7p we come to the expression 

B^u>{s) = 2eguj-^KK s [i?(s, m^) - R{ml,mjf)] + 2ieg^^KK Hs, m\). (B23) 



The imaginary part of Bp,^ is given by the tttt intermediate states so we have 

Im Bp^{s) IxnB^^is). (B24) 

Making use of (|B9|I. Bp^ takes the form 

Bpuj{s) = gp^-KTxgij^TTT, [s R(s, ml) - s R{m^p, ml) + il{s, ml)] . (B25) 

This completes our discussion of the individual functions Bij{s). The resulting inverse propagator matrix is given in 
(14. 8() . For comparison with (|3.19|l wc also list the resulting matrix elements of K 



R{ml,m,l) + -R{m%m\) 



2 Di-^Z 2 



K^^ = l-gi,^KKRi<,rni,), (B26) 
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K 

IP 


— K 




nip 


K 




= 








= bpu, - 



R{ml,ml) + -R{ml,m\) 

bpuj 9ijJ — >7rn9p — ^TTTT-^i.^p: ^tt)' 

with the functions R given in HB17|I . 

Using the values for the masses and coupHng constants from Fit III, we find the following central values for Kij-. 
Kpp = 1.017, K^^ = 0.994, Kp^ = 2.6 • lO'^, K^p = 0.2 e, K^^ = 0.06 e, e = V47rao.m. - 1/3. Clearly, the deviation 
of the matrix K from the 3x3 unit matrix is small, only at the percent level. 

Finally, we discuss the relation of our expressions to the original Gounaris-Sakurai expression for the form factor. 
The formula given in (11) of ^] can be written as 



(B27) 



where 



^GS 



fc3 



{k''[h{s) - h{ml)] + klh'{ml){ml - s)] + impVf {k/kpf nip/^s, 



(B28) 



J (i.s-m2)V2 for 
\ i{ml-\sYl^ for 



s > 4:ml, 
< s < Ami, 



(B29) 



-^log 



2K 



2m^ 



27rfc„ 



(B30) 



/is = ^log 

TT \/s 



2k 



2m, 



(B31) 



The only free parameters in the GS formula are nip and F^^. It is easy to see that the following relations hold 



ReB^p^ml) = 0, 



J-s^^BfpHs)l 



qGS 



pp 
(0) 



0, 



~C?TOpF^^. 



(B32) 



This proves (EiH) and Comparing (P28|) to l|5T2)l and ||BT7|| and setting 

1 k- 1^ i'SfW^ 



2 

.9p — 'TTTT 



247r 



1 



(B33) 



we obtain 



= 1- 



ml - s- Bpp{s) 



2 

mi 



(s-ml), 



(B34) 
(B35) 



Here Bpp{s) is defined as our Bpp{s) in HB16|I but omitting the KK contributions. 

To compare the GS with our expression it is best to choose F^{s) H4.17|l since no uj contributions are included in 
(|B35|) . Clearly, also the KK contributions arc not included in ljB35p . but one could easily do so replacing Bpp{s) by 
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the full Bpp{s). The remaining difference between (|B35|) and H4.17|l - H4.19|l is in the s-dependence of the effective p — 7 
coupling. The GS formula would correspond to the relation 

Igp^..— - 1. (B36) 

2 nip 

Using the p~ parameters from our Fit III gives, however, 

^.9p^..^ = 1.12 ±0.05. (B37) 



APPENDIX C: PROPAGATOR MATRIX AND FORM FACTORS 



In this appendix we give the derivation of the expressions H4.13|l . H4.14|l . and H4.17(l for the form factors. We start 
with writing the inverse propagator matrix (|4.8|1 as follows 







A-i(.) 



e s 



where 



and 



A-\s) = 













s bp^ + 


Bpu{s) 



s bp^ + Bp^ [s] -mf, + s + B^^ (s) 



The tranverse propagator to order e then reads 

At(s) = < 1 + e 



e s 








sA{s)G^^v 






l/s 








A(.) 



l/s 


-eG^^yAis) 


~eK{s)G^^v 


A(,s) 



+ 0(e2). 



(Cl) 



(C2) 



(C3) 



(C4) 



A(s) = (det A'\s)y 



-ml + s + B^^{s) -sbp^-Bp^{s) 



-Sbn 



BpUs) -ml + s + Bpp{s) I ' 



det A \s) = i~ml + s + B^^{s)){-ml + s + Bpp{s)) - (sbp^ + Bp^{s)f 
The pole masses and the pole widths of the p° and uj are obtained as solutions of 

det A;^\s) = s ■ det A"^(.s) = 0. 



(C5) 
(C6) 
(C7) 



To obtain the expression for the electromagnetic form factor we use now (|2.2ll . (|4.12() . (|3.11|) . and (|3.12|l . This gives 
with q = p + p' , s = (p + p')^ = 

e{p' -p)pF^{q^) = 



= {tt'^{p')'k 


-{p)\Mom 




= {■k'^{p)t:^ 


-{p)\A''mO) {-cjp.q^ + qpqu) 




~ {■k'^{p)t:^ 


-(p)||F«^0)||0)A£^)(g) 




= {Tr^{p')n^ 


-{p)\\V^^^\0)\\0)Ai^''\s) {g,pq'^ 




= {it^{p')tt' 


-ip)\\V^'lm\0)^^'^''\s) 





(C8) 
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Inserting here H4.4|l . (|4.5|l and (|C4|) leads to 

F^{s) - l~sG^^,,A{s)G,^v, (C9) 

where 

Gv^.. = ( f^"^"" ) , (CIO) 



2 



We can further simpUfy the expression for F-^ in ljC9|l taking into account that guj^-mr is much smaller than (7p_>^7rj 
and that the puj mixing parameter is also small as comes out from the fit. From ljC6|) and (|B25|) we get then 



det A-i(s) = {-ml + s + B^^{s)){-ml + s + Bpp{s)) + 0{gl^,,,bl^,g^^^^bp^) (Cll) 
For Ft^{s) of HC9|I this gives 

1 Ttj + -eB,p{s) 



TP / \ 1 I " '""^ ^ 



). (C12) 
In a similar way we get for the 77r transition form factor H2.13|l 

F^^{s) = F^,(0)-sG^_^.A(s)Gy^,^, (C13) 
with F^^(O) given by the anomaly H2.15|l and 

Expanding in gui-t-K-n and fep^^ gives 



(sj = (UJ H 2 p / X -I 2 5 — TT 

mp m-' - s - Bpp[s) - s - B^ui(s) 



sbp^ + Bp^{s) 



(m2 - s - Spp(s)) (m2 - s - B^^{s)) 



+ 0{gl^^^,bl^,g^^^^bp^). (C15) 
The weak form factor (|2.12l) F+ is obtained from ljC12l) setting gu^-mr ~ and bp^ = 0. This gives 



F^is) = - + -,9p^..^^ - + ^-^B^pis) ) . (C16) 



Using now (|B21|) leads to (|IT7jl . 



